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TECHNICAL MEMORANDUM X- 64790 


A REACTIVE TORQUE CONTROL LAW FOR GYROSCOPICALLY 
CONTROLLED SPACE VEHICLES 

SECTION I. INTRODUCTION 


The reactive torque and the well-known pseudo or generalized inverse 
methods of control represent only two of many solutions [ 1] to the same 
problem . A comparison of the two methods is helpful in introducing the 
reactive torque concept, since the two control schemes are developed by 
approaching the problem from exactly opposite viewpoints . 


The pseudoinverse is a result of the more direct approach to deriving 
a control law for gyroscopically torqued vehicles. A simplified block diagram 
of the generalized inverse scheme is shown in Figure i . The control vari- 
ables, i.e. , the gimbal rate commands to the CMG's, are determined by the 
required torques about the body-fixed axes and then making use of the fact 
that t 2] 



where 


n = 1, 2, . . . m 

T = torque applied to each of the body -fixed axes by the combined 
action of the CMG's 

<5^ = angular gimbal rate of nth CMG 



H n = momentum of the CMG rotor 

B = matrix expressing the vector cross product and resolution to 
the vehicle body axes (B is a function of fee momenta and 
gimbal angles) . 


r.) PERFORMANCE ERROR 



Figure 1. Simplified block diagram of pseudoinverse control scheme. 

Assuming the commanded (S c ) and actual (<5) gimbal rates are equal and 
solving equation (1) for the gimbal rates gives 


<5 C = B _1 T . (2) 

If there are four or more gyroscopes on board the vehicle, then B is a non- 

square matrix and many solutions for fee <5 C matrix exist. A well-known 
solution to this problem is fee pseudoinverse given by 

~ T T -1 ~ 

<5 C = B (BB ) T (3) 

This solution is usually formulated to satisfy the side condition that the 
magnitude squared of fee solution vector is a minimum . 
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The generalized inverse is a relatively complex method of control 
which requires considerable computing capability to calculate the inverse 
and to generate the torque requirements. Most of these complexities are 
circumvented by considering the control problem from the standpoint of the 
reactive torque produced on the individual CMG's by the rotating vehicle. If 
the vehicle has body rates (aJ) about the body-fixed axes, then the reactive 
torques on the CMG gimbals and opposing the torque output of the CMG servos 
are expressed by 


T = cj x H n = Mu (4) 

X\ 


where M is a matrix expressing the vector cross product and the resolution 
onto the gimbal axis of each CMG. Figure 2 shows the block diagram of the 
reactive torque control scheme. It is a classic closed-loop control system • 
operating on the error-nulling principle requiring that the outer loop be closed 
for stability. 



Figure 2. Simplified block diagram of reactive torque control scheme. 

Unlike equation (2) the reactive torques are rigorously determined by 
equation (4) , the option of satisfying some side condition does not exist, there 
is no need to generate a matrix inverse, and the rate commands (co^) are easily 
calculated. 

The reactive torque control scheme calculates the future system state 
as if the desired performance had already been achieved with the existing 
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gimbal angles. The CMG commands are then calculated on an individual basis 
and in a manner based upon the difference of the present and future states . 


SECTION II. COORDINATE SYStEMS AND CMG ORIENTATION 


Consider a vehicle controlled by four single-gimbal control moment gyro- 
scopes (SGCMG) whose momentum vectors are known with respect to the body- 
fixed axes. As a baseline configuration* , the SGCMG's are each mounted in 
the face of a pyramid with base angle 0 as shown in Figure 3. If the gimbal 
angles are defined as in Figure 3, then the positive gimbal angle rate vectors 
are in the direction shown and are perpendicular to the pyramid faces . The 
rate Vectors are measured with respect to the body axes and the CMG momen- 
tum vectors as shown in the figure represent the null position, i.e. , zero 
gimbal angle, fortheCMG's. 


x 



❖ TRW : TRW Technical Proposal Documents (Attitude Control and Determi- 
nation System) . May 1972. 
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Let pj, p 2 , p 3 , and p 4 be unit vectors which lie along the respective 
gimbal rate vectors such that 


Pj = cj3 i - sjS k 
~p<i = cj3T+ s/3 1c 

p 3 = - c/3 F+ s/3f 

p 4 = -C/3F- s/3j 


(5) 


where i, j, and k are unit vectors along the x, y, and z axes, respectively. 


If the vehicle body rates are known (from onboard sensors) and are 
defined to be 



then the vector cross product indicated in equation ( 4) becomes 

T B=“v xH n . (7) 

For future purposes, the vector multiplication as shown in equation 
(6) will not be used. The subscripts x, y, and z will be understood to denote 
scalars in the i, F> and k directions. Equation (6) will be written as 
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SECTION I IL DERIVATION OF EQUATIONS 
AND TRANSFORMATIONS 


The reactive torques have been defined as those components of co^ x 

which are parallel to the individual gimbal axes defined by the ~p vectors. 

The reactive torque vector matrix, equation (7) , is most easily computed by 
forming the vector cross products shown in Table i. 


TABLE 1. REACTIVE TORQUES BEFORE RESOLUTION 
TO CMG GIMBAL AXES 


X 

H X1 

Hx2 

Hx3 

H X4 

CO 

0 

0 

0 

0 

X 





CO 

y 

HjS^s/Sujy k 

H 2 s6 2 s/3co^Ic 

H 3 s6 3 s/3co Tc 

H 4 s6 4 si3a;^ k 

CO 

z 

-Hjsdjs^co J" 
z 

-H 2 s6 2 s/3co j* 
z 

-H 3 s6 3 s/Ico f 
z 

-H 4 s6 4 s/?co J 
z 



— 


: - • 


M 

Hyi 

Hy 2 

% 3 

H y* 

CO 

X 

-ILcdiCo k 
1 1 X 

HocdoOJ k 
* * X 

-H 3 s<5 3 cj3co k 

X 

H 4 s<5 4 cj3co x k 

CO 

0 

0 

0 

0 

y 





CO 

z 

ELcSiCO i 
1 1 z 

-H 2 c<5 2 co i 
z 

HgSSgcPc^r 

-H 4 S<5 4 c|3co i 
z 


X 

Hzi 

. H Z2 

Hz 3 

H Z4 

CO 

X 

H^sSjC/Sco J" 

X 

-H 2 s5 2 c/3co f 

X 

H 3 c6 3 w x f 

-H 4 c5 4 co x f 

CO 

-HiSdiCjSco F 

H 2 sa9Ci3co r 

-H^cdfiCO F 

H 4 c6 4 co i 

y 

y 

" ‘ y 

y 

y 

CO 

z 

0 

0 

0 

0 


Note: The table is read in the following manner. The vector cross product of 
co and H z4 is co x H Z4 = H 4 c6 4 co T 

y y y 
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The values of the table are calculated by writing the CMG momenta in 
vector notation as 


H - H i+H j + H Tc 
n xn yn J zn 


(9) 


and then forming the cross product of o> v with each of the momentum compo- 
nents of each of the four SGCMG's. The values shown in the table represent 
all reactive torques before resolution to the CMG gimbal axes. 

Using Table 1 and equation (5) , the reactive torques on the individual 
SGCMG's produced by the rotating vehicle and opposing the torques driving the 
gimbals can be expressed by 


h E1 \ A<“ v xH i) • PiipA 


T = 
R 


T 


R2 


R3 


l I 

\ 7 / 

\ R4/ 


Ku x H 2 ) • p 2 ]p 2 


[(w y xH 3 ) • p 3 ]p 3 


(K) ( w y ) 


[(co v X h 4 ) • p 4 ]p 4 


/ 


( 10 ) 


where K is the matrix defined by 


(K) = 


^HiCSiS/3 

H 2 c6 2 s/3 
H 3 cS 3 s/3 
H 4 c5 4 s/3 




-HjSS! 

H 2 s5 2 
H 3 c6 3 c p 
-H 4 c6 4 c/3 



(ID 


The reactive torque vectors as defined by equation ( 3) can best be 
visualized by examining the vector diagram of an individual SGCMG. Figure 
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4 represents the reactive torques and body rates as seen by SGCMG number i 
when o < < 71/2 radians. 


x 



" ' Y 

Figure 4. Reactive torques and body rates as seen by CMG number i. 
From Table i 



Since T ly is perpendicular to p t , it contributes nothing to the value of . 
From Figure 4 the magnitude of T is 
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( 13 ) 


which is identical to 

T^V^-Pl ‘ < 14 > 

SECTION IV. SYNTHESIS OF THE STEERING LAW 

Assume the driving torques on the gimbals to be produced by rate- 
commanded servos whose outputs are the torques defined by 



Since the CMG's are single gimbal, these torques are fixed in direction 
relative to the body axes and lie along the gimbal axes defined by the unit p 
vectors . 

Let the gimbal inertias (which may or may not be identical) be defined 
as that inertia "seen” by the rate commanded servos. Symbolically the 
gimbal inertias will be defined by the equation 




then 


* 

If <5 now represents the angular gimbal rates as shown in Figure 3, 
n 



where S is the Laplace operator and, as used, signifies the differentiation of 
8 with respect to time. 

Equation (18) offers a convenient method of commanding the vehicle 
to rotate in some prescribed manner. It is reasonable to assume that the 
commanded and the actual gimbal rates are essentially equal. If the 
commanded gimbal rates are 
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then equation (18) becomes 


6 = (si r 1 

c g 


(T m-V 


( 20 ) 


Equation 20 defines the manner in which the torque requirements will 

be distributed among the various CMG's. If T is allowed to change by some 

ft 

amount A T R then the control law forces a corresponding change in the gimbal 
rate command according to equation ( 20) . 

Let 

VVat r ’ (21) 

then substitution into equation (20) gives 




8 + A 5 = (SI ) -I (T - T - A T ) 
c c g m R R 


- (SI r 1 (T - T ) - (SI r 1 (AT) 
g m R g R 


It now follows that 


A 8 = - (S I ) -1 (AT) 
c g R 


( 22 ) 


Assume how that it is desired to command a change in the body rate 

vector co . Let the commanded body rates be 
v 


/ “ xg \ 


(co ) = 


CO 


yg 


(23) 


\ 


CO 


Zg / 
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For a given change in the vehicle body rates the corresponding change 
in the reactive torques can be calculated from equation (10) as follows 

AT = (K)(w ) - (K)(w ) = (K) (« -Z) . (24) 

R g v g v 

Let the rate error matrix be 

(«,) = (w y - w g ) . (25) 

The substitution of equations (24) and (25) into (22) gives the result 

AS = (SI ) -1 (K) (w ) . (26) 

eg e 


Equation (26) calculates the necessary adjustment in the gimbal rate commands 

to cause the body rate errors (co ) to approach zero. 

6 

No provision has been made at this point to determine the orientation 
of the vehicle. For this reason, a more general formulation of equation (26) 
is 

A 6 = (SI )~ i (K) (co ) (27) 

e g c 


If the initial values (6 .) of the commanded gimbal rates are known then 

Cl 

y =y. + A6 r . 

C Cl c 


The purpose of the introduction of the u> c matrix is to make the system 
dependent upon the attitude error as well as the attitude ra.te error defined by 
equation ( 25) . Since small signal orientation errors of the vehicle can be 


12 



determined by integrating the body rates, the position error matrix (e ) can 
also be defined in terms of equation (25) f ' 



If the ^attitude rate and position control gains are also included in the 

definition of co then 
c 


“ = (ai) (w ) + (a ) (0 ) = 
c e o e 



(29) 


where 



The control equation (27) used in conjunction with equation (29) pro- 
duces the simplified system block diagram shown in Figure 5. The individual 
blocks have been divided into several major areas by the dashed lines for the 
purpose of a clearer definition or to indicate where the functions are performed 
in the spacecraft. 
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Figure 5 . Simplified block diagram of reactive torque control scheme . 
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SECTION V, COMMENTS AND CONCLUSIONS 

Single-gimbal control moment gyroscopes were selected for this report 
in order to keep the equations as simple as possible and still illustrate the 
control principles involved . The reactive torque-control method will also 
apply to double gimbal CMG's and for any reasonable arrangement of the 
individual gyroscopes. 

While it is doubtful that the gyroscopes are commanded in a way which 
minimizes the sum of the individual gimbal rates (as in the pseudoinverse 
method) , the reactive torque method does command the gimbals in a gradient 
manner which gives the smallest commands to the least effective gyroscopes . 

Normally, the vehicle body rates will be determined by rate gyroscopes 
In order to null any errors which may develop, onboard updates of attitude 
rate and position may be accomplished by means of star trackers and sun 
sensors . The outputs of these sensors would be used to modify the rate gyro 
outputs and to determine the orientation of the vehicle . 

Singularity avoidance has not been discussed in this report. Singular- 
ities do exist and should be the next point of investigation for anyone consider- 
ing implementing this method, of control. 

The reactive torque control law has been simulated only to the point of 
showing that it is a feasible method of controlling a space vehicle. It repre- 
sents a very simple system when compared to most existing methods and 
should be easily implemented . 
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